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ABSTRACT 

We present photometry with the Advanced Camera for Surveys (ACS) on 
the Hubble Space Telescope (HST) of stars in the Magellanic starburst galaxy 
NGC 4449. The galaxy has been imaged in the F435W (B), F555W (V) and 
F814W (I) broad-band filters, and in the F658N (Ha) narrow-band filter. Our 
photometry includes ~ 300,000 objects in the (B, V) color-magnitude diagram 
(CMD) down to V < 28, and ^ 400,000 objects in the (V, I) CMD, down to I 
< 27 . A subsample of ~ 200,000 stars has been photometrized in all the three 
bands simultaneously. The features observed in the CMDs imply a variety of 
stellar ages up to at least 1 Gyr, and possibly as old as a Hubble time. The spatial 
variation of the CMD morphology and of the red giant branch colors point toward 
the presence of an age gradient: young and intermediate-age stars tend to be 
concentrated toward the galactic center, while old stars are present everywhere. 
The spatial variation in the average luminosity of carbon stars suggests that there 
is not a strong metallicity gradient (< 0.2 dex). Also, we detect an interesting 
resolved star cluster on the West side of the galaxy, surrounded by a symmetric 
tidal or spiral feature consisting of young stars. The positions of the stars in 
NGC 4449 younger than 10 Myr are strongly correlated with the Ha emission. 
We derive the distance of NGC 4449 from the tip of the red giant branch to 
be D = 3.82 ± 0.27 Mpc. This result is in agreement with the distance that we 
derive from the luminosity of the carbon stars. 
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Introduction 



Starbursts are short and intense episodes of star formation (SF) that usually occur in the 
central regions of galaxies and dominate their integrated light. The associated star- formation 
rates (SFR) are so high that the existing gas supply can sustain the stellar production only 
on timescales much shorter than a cosmic time (< 1 Gyr). 

The importance of the starburst phenomenon in the context of cosmology and galaxy 
evolution has been dramatically boosted in recent years by deep imaging and spectroscopic 
surveys which have discovered star-forming galaxies at high redshift: a population of dusty 
and massive starbursts, with SFRs as high as ~ 100 - 1000 Mp^ y^|'~^i has been unve i led in 
the submillimeter and millimeter wavelengths at z>2 (IBlain et al.ll2002l : IScott et al.ll2002l ) 
and star-fo rming galaxies at z > 3 have been di scovered with the Lyman break selection 



technique (ISteidel et al.l Il996l : iPettini et al.l 120011) and thr ough Lyman-a emission surveys 



fle.g.. Rhoads et all J200of )7 see also iLe Fevre et al.l (120051 ) for a more recent independent 



approach) . 

In the local Universe, starbursts are mostly found in dwarf irregular galaxies, and con- 
tribute ~ 2 5% of the whole massive SF (|Heckmarull998r). Both observations and theoreti- 
cal models (ILarson fc Tinsleylll978l : le.g.. Genzel et al.lll998l : le.g.. Noguchi fc Ishibashilll986l ) 
show that strong starbursts are usually triggered by processes such as interaction or merging 
of galaxies, or by accretion of gas, which probably played an important role in the formation 
and evolution of galaxies at high redshift. Thus, nearby starbursts can serve as local analogs 
to primeval galaxies to test our ideas about SF, evolution of massive stars, and physics of the 
interstellar medium (ISM) in "extreme" environments. The high spatial resolution and high 
sensitivity of Hubble Space Telescope offer the possibility to study the evolution of nearby 
starbursts in details. This is fundamental in order to address many of the still open questions 
in cosmological astrophysics: What are the main characteristics of primeval galaxies? What 
is the nature of star-forming galaxies at high redshift? How important are accretion and 
merging processes in the formation and evolution of galaxies? 

The Magellanic irregular galaxy NGC 4449 (asooo = 12'^28'"11^9, ^2000 = +44°05'40", 
/ = 136.84 and b = 72.4), at a distance of 3.82 ± 0.27 Mpc (see Section 5), is one of 
the best studied and spectacular nearby starbursts. It has been observed across the whole 
electromagnetic spectrum and displays both interesting and uncommon properties. It is one 
of the most luminous and active irregular galaxies. Its integrated magnitud e Mr = — 18. 2 
makes it ^ 1.4 times as luminous as the Large Magellanic Cloud (LMC) (IHunteii 119971 ). 
Thronson et al.l (119871 ) estimated a current SFR of ~ 1.5 M© yr^^. NGC 4449 is also 



the only local example of a global st arburst, in the sense that the current SF is occurring 
throughout the galaxy (IHunteii 119971 ). This makes NGC 4449 more similar to Lyman break 
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Galaxies (LBGs) at high redshift (2; ^ 3), where the brightest regions of SF are embedded 
i n a more diffuse nebulosity and dominate the integrated light also at optical wavelengths 
jGiavaliscol 120021 ) ■ 



Abundance est imate s in NGC 444 9 were derived in the HII regio ns bv iTalent 



Hunter et all (Il982f ) and [Martini (Il997h . and for NGC 4449 nucleus bv iBoker et al 



1980 ) 



(120011 ) ■ 



The published values are in good agreement wit h each other, and provide 12 + log(0/H) ^ 
8.31. Adopting the oxygen solar abundance from lGrevesse fc Sauvall (119981 ). 12 + log(0/H)Q = 
8.83, we obtain [0/H] = -0.52, i.e. NGC 4449 oxygen content is almost one third solar, as in 
the LMC. New solar abundance estimates, based on 3D hydrodynamic models of the solar 
atmosphere, accounting for departures from LTE, arid on improved atomic and molecular 
data, provide 12 + log(O/H)0 = 8.66 (lAsplund et al.ll2005l ). However, the new lower abun- 
dances seem to be inconsistent with helioseismo logy data, u nless the majority of the inputs 
needed to make the solar model ar e chan ged (IBasul 120071 ). Thus, we will adopt the old 
abundances from iGrevesse fc Sauvall (119981 ) throughout the paper. 



Radio observations of NGC 4449 have shown a very extended HI halo (~ 90 kpc in 
diameter) which is a factor of ~ 10 larger than the optical diamete r of the galaxy an d 
appears to rotate in the opposite direction to the gas in the center (jBajaja et al.l Il994j ). 
Hunter et al. (1998, 1999) have resolved this halo into a central disk-like feature and large 
gas streamers that wrap around the galaxy. Both the morphology and the dynamics of 
the HI gas suggest that NGC 4449 has undergone some interaction in the past. A gas-rich 
companion galaxy, DD O 125, at the projected distance of ~ 40 kpc, could have been involved 
kheis fc Kohl]l200ll ). 



NGC 4449 has numerous (~ 60) star clusters with ages up to 1 Gyr (iGelatt et al.ll200ll ) 
and a young (~ 6-10 Myr) central cluster (iBoker et al.ll2001l). a prominent stellar bar which 
covers a large fraction of the optical b ody (IHunter et al 



19991 ). and a spherical distribution 



of older (3-5 Gyr) stars (|Bothunlll986l ). The galaxy has also been demonstrated to contain 
molecular clouds from CO observations (iHunter &: Thronsonlll996f) and to have an infrared 



10-150 /im) luminosity of 2 x 10 erg s ( iThronson et al.lll987l ). The ionized gas shows a 



very turbulent morphology with filaments, shells and bubbles which extend for several kpc 
(Hunter & Gallagher 1990, 1997). The kinematics of the HII regions within the galaxy is 



chaotic, again suggesting the possibility of a collision or merger (le.g.. Valdez-Gutierrez et al. 



2OO2I ). Some 40% of the X-ray emission in NGC 4449 comes from hot gas with a complex 
morphology similar to t hat observed in Ha, implying an expanding super-bubble with a 
velocity of ~ 220 kms~^ ( ISummers et al.ll2003l ). 



All these observational data suggest that the late-type galaxy NGC 4449 may be chang- 
ing as a result of an external perturbation, i.e., interaction or merger with another galaxy. 
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or accretion of a gas cloud. A detailed study of the star-formation history (SFH) of this 
galaxy is fundamental in order to derive a coherent picture for its evolution, and understand 
the connection between possible merging/accretion processes and the global starburst. With 
the aim of inferring its SFH, we have observed NGC 4449 with the Advanced Camera for 
Surveys (ACS) on the Hubble Space Telescope (HST) in the F435W, F555W, F814W and 
F658N filters. In this paper we present the new data and the resulting color-magnitude 
diagrams (CMDs) (Sections 2, 3 and 4). We derive a new estimate of the distance modulus 
from the magnitude of the tip of the red giant branch (TRGB) and the average magnitude 
of the carbon stars in Section 5. In Section 6, the empirical CMDs are compared with stellar 
evolutionary tracks. With the use of the tracks, we are able to derive the spatial distribution 
of stars of different age in the field of NGC 4449. The Conclusions are presented in Section 7. 
The detailed SFH of NGC 4449 will be derived through synthetic CMDs in a forthcoming 
paper. 



2. Observations and data reduction 

The observations were performed in November 2005 with the ACS Wide Field Camera 
(WFC) using the F435W (B), F555W (V) and F814W (I) broad-band fihers, and the F658N 
(Ha) narrow-band filter (GO program 10585, PI Aloisi). We had two different pointings in a 
rectangular shape along the major axis of the galaxy. Each pointing was organized with a 4 
- exposure half -|- integer pixel dither pattern with the following offsets in arcseconds: (0,0) 
for exposure 1, (0.12, 0.08) for exposure 2, (0.25, 2.98) for exposure 3, and (0.37, 3.07) for 
exposure 4. This dither pattern is suitable to remove cosmic rays and hot /bad pixels, fill the 
gap between the two CCDs of the ACS/WFC, and improve the PSF sampling. The exposure 
times in the different broad-band filters were chosen to reach at least one magnitude fainter 
than the TRGB, at I f»24, with a photometric error below 0.1 mag (S/N > 10). Eight 
exposures of ~ 900 s, 600 s, 500 s and 90 s were acquired for each of the B, V, I and Ha 
filters, respectively. 

For each filter, the eight dithered frames, calibrated through the most up-to-date version 
of the ACS calibration pipeline (CALACS), w ere co- added into a sin gle mosaicked image 



using the software package MULTIDRIZZLE (IKoekemoer et al.ll2002l ). During the image 
combination, we fine-tuned the image alignment, accounting for shifts, rotations, and scale 
variations between images. The MULTIDRIZZLE procedure also corrects the ACS images 
for geometric distortion and provides removal of cosmic rays and bad pixels. The total 
field of view of the resampled mosaicked image is ~ 380 x 200 arcsec^, with a pixel size of 
0.035 "(0.7 times the original ACS/WFC pixel size). 
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To choose the optimal drizzle parameters, we experimented with different combinations 
of the MULTIDRIZZLE parameters pixfrac (the linear size of the "drop" in the input pixels) 
and pixscale (the size of output pixels). One must choose a pixfrac value that is small enough 
to avoid degrading the final image, but large enough that, when all images are dropped on 
the final frame, the flux coverage of the output image is fairly uniform. Statistics performed 
on the final drizzled weight image should yield an rms value which is less than 20% of the 
median value. In general, the pixfrac should be slightly larger than the scale value to allow 
some of the 'drop' to spill over to adjacent pixels. Following these guidelines, we find that 
pixfrac=0.8 and pixscale=0.7 provide the best resolution and PSF sampling for our dithered 
images. 

The total integration times are ~ 3600 s, 2400 s, 2000 s and 360 s for the B, V, I and 
Ha images, respectively. Only in a small region of overlap between the two pointings (~ 
30 X 200 arcsec^) the integration times are twice as those listed above. Figure [1] shows the 
mosaicked true-color image created by combining the data in the four filters. 

The phot ometric reduction of the images was performed with the DAOPHOT package 
( Stetson 1987 ) in the IRAF environmenll^. The instrumental magnitudes were estimated via 



a PSF-fitting technique. We constructed a PSF template for each of the four ACS chips 
(2 CCDs at 2 pointings) contributing to the final mosaicked image. To derive the PSF, 
we selected the most isolated and clean stars, uniformly distributed within each chip. The 
PSF is modeled with an analytic moffat function plus additive corrections derived from the 
residuals of the fit to the PSF stars. The additive corrections include the first and second 
order derivatives of the PSF with respect to the X and Y positions in the image. This 
procedure allows us to properly model the spatial variation of the PSF in the ACS/WFC 
field of view. 

The stars were detected independently in the three bands, without forcing in the shal- 
lowest frames the detection of the objects found in the deepest one. For comparison, we also 
ran the photometry on a list of stars detected on the sum of the B, V and I images. The 
luminosity functions (LFs) obtained in the various bands with the two different approaches 
are presented in Fig. [2J We notice that the forced search pushes the detection of stars ~ 0.5 
magnitude deeper than the independent search. On the other hand, upon closer inspection 
the majority of the "gained" objects turn out to be spurious detections or stars with large 
photometric errors. Furthermore, the deeper photometry is not deep enough to detect the 
next features of interest in the CMD (the horizontal branch, the red clump or the asymp- 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by AURA, 
Inc., under cooperative agreement with the National Science Foundation 
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totic giant branch (AGB) bump, see Section 3), and thus it does not provide any additional 
information for our study. In the following we thus use the photometry obtained with the 
independent search on the three bands. 

Aperture photometry with PHOT, and then PSF-fitting photometry with the ALLSTAR 
package, were performed at the position of the objects detected in the B, V and I images. 
The instrumental magnitudes were measured adopting the appropriate PSF model to fit the 
stars according to their position in the frame. The B,V and I catalogs were then cross- 
correlated with the requirement of a spatial offset smaller than 1 pixel between the positions 
of the stars in the different frames. This led to 299,115 objects having a measured magnitude 
in both B and V, 402,136 objects in V and I, and 213,187 objects photometrized in all the 
three bands simultaneously. 

The conversion of the instrumental mag nitudes m,; to t h e HST VEGAMAG system was 



performed by following the prescriptions in ISirianni et al.l (120051 ). The HST VEGAMAG 



magnitudes are derived according to the equation: 



m = mi + Cap + Cinf + ZP + CcTE, (1) 

where is the DAOPHOT magnitude (—2.5 x log(counts/exptime)) within a circular aper- 
ture of 2-pixel radius, and with the sky value computed in an annulus from 8 to 10 pixels; 
Cap is the aperture correction to convert the photometry from the 2 pixel to the conventional 
0.5" radius aperture, and with the sky computed at "infinite" ; we compu ted Cap from iso- 



lated stars selected in our images; Cinf is taken from ISirianni et al.l (120051 ) and is an offset 
to convert the magnitude from the 0.5 "radius into a nominal infinite aperture; ZP is the 
HST VEGAMAG zeropoint for the given filter. Corrections for imperfect charge transfer 
efficiency (GTE) were calculated from each single expousure, and then averaged, following 



the formulation of iRiess fc Mackl (120041 ) (Eq. [2]), which accounts for the time dependence of 



the photometric losses: 

CcTE = 10^ X SKY^ X FLf/X^ x ^ x ^^^-^ " 52333 

2048 365 ^ ^ 

where SKY is the sky counts per pixel per exposure, FLUX is the star counts per expo- 
sure within our adopted photometry aperture (r=2 pixel in the resampled drizzled image, 
corresponding to 1.4 pixel in the original scale), Y is the number of charge transfers, and 
MJD is the Modified Julian Date. The coefficients of equation ([2]) were extrapolated from 
a r=3 pixel aperture to a r=1.4 pixel aperture (in the original scale), and are A = 1.08, 
B = —0.309, and C = —0.976. The computed GTE corrections are negligible for the bright- 
est stars, but can be as high as ~ 0.1 mag for the faintest stars. We did not transform the 
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final magnitudes to the Johnson- Cousins B, V, I system, since this would introduce addi- 
tional uncertainties. How ever, such transform ations can be done in straightforward manner 



using the prescriptions of ISirianni et al.l (120051 ) . The ACS VEGA MAG magnitudes we re not 



corrected for Galactic foreground extinction (E(B — V) = 0.019, ( Schlegel et al. IQQsl )) and 



internal reddening. Concerning the internal reddening. Hill et al. ( IQQsl ) derived E(B — V) 



0.18 from the ifa/H/3 ratio measured in NGC 4449 HII regions. This value can be considered 
an upper limit to the average internal extinction, since the nebular gas is usually associated 
with young star forming regions, which tend to be inherently more dusty than th e regions 



i n wh ich older stars reside (as demonstrated explicitly for the case of the LMC; IZaritsky 



f|l999h l. 



Because of both the large number of dithered exposures and the conservative approach 
of the independent search on the three images, our catalog is essentially free of instrumental 
artifacts such as cosmic rays or hot pixels. The distribution of the DAOPHOT parameters 
cr, and sharpness is shown in Fig. [3l The a parameter measures the uncertainty on the 
magnitude, the is the residual per degree of freedom of the PSF-fitting procedure, and the 
s/iarpnes5 provides a measure of the intrinsic size of the object with respect to the PSF. Notice 
that the DAOPHOT/ ALLSTAR package automatically rejects the objects with a > 0.55. 
The obtained distributions for a, and sharpness suggest that the vast majority of the 
detected sources are stars in the galaxy with a small contamination from stellar blends and 
background galaxies. There are some objects with very bright magnitudes (mF8i4w ^ 22) 
and sharpness >0.5 in Fig. [3] (this is also observed in the F435W and F555W filters, for which 
we do not show the and sharpness distributions). Their sharpness values imply that they 
have a larger intrinsic size than the PSF, and thus may not be individual stars. By visually 
inspecting these objects in all the images, we recognized several candidate star clusters and 
background galaxies. Some of the candidate star clusters look like fairly round but extended 
objects; some others present a central core, and are partially resolved into individual stars 
in the outskirts. We detect at least 42 clusters in our data, some of which look like very 
massive globular clusters. The candidate clusters and the galaxies that were identified by 
eye were rejected from the photometric catalog. We are left with 299,014 objects in the (B, 
V) catalog, 402,045 objects in the (V, I) catalog, and 213,099 objects photometrized in all 
the three bands. We experimented with many other cuts in a, and sharpness, but none 
of them affected the global appearance of the CMDs and the detected evolutionary features. 
A detailed study of the cluster properties will be presented in a forthcoming paper (Aloisi 
et al., in preparation). 
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3. Incompleteness and blending 

To evaluate the role of incompleteness and blending in our data, we performed arti- 
fici al star experimen ts on the drizzle-combined frames, following the procedure described 



by iTosi et al.l (1200 ll ). These tests serve to probe observational effects associated with the 
data reduction process, such as the accuracy of the photometric measurements, the crowd- 
ing conditions, and the ability of the PSF-fitting procedure in resolving partially overlapped 
sources. We performed the tests using to the following procedure. We divided the frames 
into grids of cells of chosen width (50 pixels) and randomly added one artificial star per 
cell at each run. This procedure prevents the artificial stars to interfere with each other, 
and avoids to bias the experiments towards an artificial crowding not really present in the 
original frames. The position of the grid is randomly changed at each run, and after a large 
number of experiments the stars are uniformly distributed over the frame. In each filter, we 
assign to the artificial star a random input magnitude between mi and m2, with mi ~ 3 mag 
brighter than the brightest star in the CMD, and m2 ~ 3 mag fainter then the faintest star in 
the CMD. At each run, the frame is re-reduced following exactly the same procedure as for 
the real data. The output photometric catalog is cross-correlated with a sum of the original 
photometric catalog of real stars and the list of the artificial stars added into the frame. 
This prevents cross-correlation of artificial stars in the input list with real stars recovered in 
the output photometric catalog. We simulated about half a million stars for each filter. At 
each magnitude level, the completeness of our photometry is computed as the ratio of the 
number of recovered artificial stars over the number of added ones. The completeness levels 
in the color magnitude diagrams (see Section 4) are the product of the completeness factors 
in the two involved passbands. 

We show in Fig. H] the Am difference between the input and output magnitudes of the 
artificial stars as a function of the input magnitude, for the F435W, F555W and F814W fil- 
ters. The solid lines superimposed on the artificial star distributions correspond to the mean 
Am (central line), and the ±1(7^ values around the mean. The plotted Am distributions 
provide a complete and statistically robust characterization of the photometric error as a 
function of magnitude, for each filter. By comparing the am with the DAOPHOT errors in 
Fig. [3|, it is apparent that the DAOPHOT package increasingly underestimates the actual 
errors toward fainter magnitudes. For instance, for a star with V~ 25.5 and I~ 24, (tip 
of the red giant branch, see Section 5), the mean DAOPHOT error is ~ 0.05 mag in both 
bands, while the from the artificial star tests is ~ 0.15 and ~ 0.1 in V and I, respectively. 
The systematic deviation from of the mean Am indicates the increasing effect of blending, 
i.e. faint artificial stars recovered brighter than in input because they happen to overlap 
other faint objects. 
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4. Color-magnitude diagrams 

The CMDs are shown in Figures [5] and [61 We plot the nipsssw versus mF435w — mpsssw 
CMD of the 299,014 stars matched between the B and V catalogs in Fig. [5|, and we plot the 
mF8i4W versus nipsssw ~ mF8i4W CMD of the 402,045 stars matched between the V and I 
catalogs in Fig. O We indicate the 90 % (solid line) and 50 % (dashed line) completeness 
levels as derived from the artificial star experiments on the two CMDs. The average size 
of the photometric errors at different magnitudes, as derived from artificial star tests, is 
indicated as well. 

The two CMDs show all the evolutionary features expected at the magnitudes sampled 
by our data: a well defined blue plume and red plume, the red horizontal tail of the carbon 
stars in the mF8i4w versus mpsssw — niF8i4w CMD, and a prominent red giant branch (RGB). 
The blue plume is located at mF435w — ^F555W and mpsssw — niF8i4w — —0.1 in the two 
diagrams, with the brightest stars detected at mpsssw? mF8i4W ~ 18- It samples both stars 
in the main-sequence (MS) evolutionary phase and evolved stars at the hot edge of the core 
helium burning phase. The blue plume extends down to the faintest magnitudes in our 
data, at mF555w ~ 28. The red plume is slightly inclined with respect to the blue plume, 
with mF555w ^ 25 and colors extending from mF435w — niF555w ~ 1.4 to ~ 1.8 in the 
inF555w, niF435w " niF555w CMD, and mpsMW ^ 23.5 and colors extending from mF555w — 
mF8i4w ~ 1-4 to ~ 2.2 in the mF8i4w versus mF555w — mF8i4w CMD. It is populated by red 
supergiants (RSGs) at the brighter magnitudes, and AGB stars at fainter luminosities. At 
intermediate colors, below mFsssw ~ 25 and mF8i4w ~ 23.5, we recognize the blue loops of 
intermediate-mass stars in the core helium burning phase. The concentration of red stars at 
K1F555W ^ 25.5 and mF8i4w ^ 24, corresponds to low-mass old stars in the RGB evolutionary 
phase. Finally, a pronounced horizontal feature, at mF8i4w ~ 23.5, and with colors extending 
from mF555w — niF8i4w ~ 1.8 to as much as mF555w — niF8i4w ~ 4, is observed in the mF8i4w, 
niF555w — niF8i4w CMD. This red tail is produced by carbon stars in the thermally pulsing 
asymptotic giant branch (TP-AGB) phase. 

In order to reveal spatial differences in the stellar population of NGC 4449, we have 
divided the galaxy's field of view into 28 (7 x 4) rectangular regions, as shown in Fig. [71 
The size of the regions (~ 55 x 55 arcsec^, corresponding to ^ 1 x 1 kpc^ at the distance 
of NGC 4449) allows us to follow spatial variations at the kpc scale, being at the same 
time large enough to provide a good sampling. The mpsssw? mF435W — i^fsssw; and mF8i4W! 
mF555w — niF8i4w CMDs derived for the different regions are shown in Fig. [HI and [H respec- 
tively. The completeness levels plotted on the CMDs of the central column show that the 
photometry is deeper in the external regions than in the galaxy center, where the high crowd- 
ing level makes the detection of faint objects more difficult. The errors, as estimated from 
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the artificial star experiments, increase toward the galaxy center, as an effect of the higher 
crowding level and the higher background. Bright stars are mostly concentrated toward 
the galaxy center, and only a few of them are present at large galactocentric di stances, in 



agreem ent with what was already observed in other dwarf irregular galaxies, (e.g.. iTosi et al. 
( l200ll )). The external region (6,4) makes the exception to this observed global trend, show- 
ing a prominent blue plume in both the CMDs. The luminous blue stars observed in these 
CMDs correspond in the image of Fig. [1] to a symmetric structure. This structure is more 
clearly visible in the top right of Fig. [111 which will be discussed in Section 6 below. This 
structure could be due to tidal tails or spiral-like feature associated with a dwarf galaxy that 
is currently being disrupted. The structure is centered on a resolved cluster-like object that 
could be the remnant nucleus of this galaxy. Fig. [7] shows a blow-up of this object. 



4.1. Carbon stars 



In the mF8i4W; nipsssw — niF8i4w CMD of Fig. El the horizontal red tail, at magnitudes 
brighter than the TRGB, is due to carbon-rich stars in the TP-AGB phase. Since AGB 
stars trace the stellar populations from~ 0.1 to several Gyrs , this well defined feature is 
suitable to investigate the SFH from old to intermediate ages ( Cioni et al. I2OO6 ). Rece ntly, 
theoretical models of TP-AGB stars have been presented by iMarigo &: Girardil (120071 ) for 
initial masses between 0.5 and 5.0 Mq and for different metallicities. Their Fig. 20 shows 
that the position of the carbon-star tracks in the log L/L© vs logTe// plane depends on both 
mass and metallicity. Higher-mass stars exhibit larger luminosities for a given metallicity. 
On the other hand, lower metallicities imply both a wider range of masses undergoing the 
C-rich phase, and higher luminosities for the more massive carbon stars. 

From an empirical aspect, we c an get the dependence of th e carbon-star luminosity on 
age and metallicity from the work of iBattinelli &: Demerd (120051 ) (hereafter BD05). The au- 
thors provide the following relation for the dependence of the carbon-star I band magnitude 
on metallicity: 



< Mi,earbon >= "4.33 + 0.28 X [Fc/H], (3) 

which was derived through a least-square fit of the mean absolute 1 band magnitude of 
carbon stars in nearby galaxies with metallicities —2 < [Fe/H] < —0.5. According to ([3|), a 
drop of 1 dex in metallicity results in a decrease of ~ 0.3 in the average magnitude of carbon 
stars. The age dependence of the carbon-star luminosity is more difficult to quantify since it 
requires a detailed knowledge of the SFH in galaxies. BD05 do not study such a dependence, 
but some qualitative considerations can be obtained from examination of their Fig. 4. The 
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scatter of the < Mi^carbon > versus [Fe/H] relation is of the order of 0.1 mag (excluding Andll 
and AndVII from the fit), i.e. ~ 20 % of the variation in < Mi carbon > spanned by the data. 
We also notice that at fixed metallicity, galaxies with current star formation (empty dots in 
Fig. 1 of BD05) tend to have brighter < Mi carbon > values, while galaxies with no current 
star formation preferentially lie below the relation. This suggests that at least part of the 
scatter of relation ([3]) is due to a spread in age, with younger stellar populations having 
bri ghter carbon stars t han o lder stellar populations. This age dependence is consistent with 



the iMarigo fc Girardil (120071 ) models, where more massive (younger) carbon stars tend to be 



more luminous. 

Carbon stars were selected in NGC 4449 at 23 < mF8i4w < 24, msssw — niF8i4w > 2.4, 
for each of the 28 regions shown in Fig. [71 The color limit was chosen to avoid a significant 
contribution of RGB and oxygen-rich AGE stars. For each region, we fitted the mF8i4w- 
band LF with a Gaussian curve, and adopted the peak of the best-fitting Gaussian as the 
average < mi carbon > in that region. The errors on the results of the Gaussian fits can 
be approximated as Am ^ a/^/N, where a is the width of the Gaussian and N is the 
number of stars. We experimented with different cuts in color (up to msssw — niF8i4w > 2.8) 
and different binnings of the data, and found that none of them significantly affects the 
final results of the analysis. The results presented in Fig. [10] were obtained by binning 
the stellar magnitudes in bins of 0.2 mag. The quantity Ami c along the ordinate is the 
difference between the carbon-star magnitude measured in a specific region, and the carbon- 
star magnitude averaged over the whole field of view of NGC 4449. Along the abscissa is 
the X coordinate in pixels. From top to bottom, the panels refer to regions of decreasing Y 
coordinate. Fig. [10] shows that the observed carbon-star luminosity is reasonably constant 
over almost the whole galaxy, within the errors. The only significant variation is observed 
in the central regions (3000 < X < 8000, panel c), where the carbon stars appear to be up 
to ~ 0.25 mag brighter than the average value. 

We performed Monte Carlo simulations to understand if the observed variation of the C 
star LF is intrinsic, i.e. due to stellar population gradients present in NGC 4449 field, or if 
it is an effect of the photometric error and completeness level variations across the field. We 
adopted an estimated intrinsic LF for the C stars, that we assumed to be constant all over 
the field. Then we investigated how the intrinsic LF is transformed into the observed one 
after completeness and photometric errors were applied at different galacto-centric distances. 
For simplicity, we assumed that the intrinsic mF8i4 distribution of the C stars is Gaussian, 
with parameters (mF8i4w,o = 23.64, a ~ 0.3) derived from the Gaussian fit to the observed 
LF in the most external regions ((1,1:4), (1:7,4), (7,1:4); where this notation will indicate the 
some of these regions). We also assumed that the C stars are not homogeneously distributed 
in color, but follow a power law in mFsssw ~ niF8i4w- The power law's parameters were 



derived by fitting the observed C star color distribution in the same external regions. Our 
assumption that the observed outer LF is a good description of the intrinsic one in those 
regions is reasonable, since the most external mF8i4, nipsssw ~ mF8i4W CMDs are more than 
90 % complete at the average C star luminosity, and the errors apsuw are only a few hundreds 
of magnitudes there (see Fig. [9]). Monte Carlo extractions of (mF8i4, mpsssw — niF8i4w) pairs 
were drawn from the assumed magnitude and color distributions. A completeness and a 
photometric error were then applied to each extracted star, for both an external and an 
internal region of NGC 4449. The effect of incompleteness and photometric errors on the 
distribution is shown in the bottom panel of Fig. [11] In the top panel we show instead the 
observed LFs for an external region ((1,1:4), (1:7,4), (7,1:4)), and for the most internal one 
(4,3), which displays the largest shift of the distribution peak with respect to the external 
region. Our simulations in the bottom panel show that the C star LF is mostly unaffected 
by incompleteness and photometric errors in the most external regions. This result testifies 
that the observed LF is very close to the intrinsic one in the periphery, and that we adopted 
a reasonable input distribution for the Monte Carlo simulations. In the most internal region, 
instead, the LF is shifted toward brighter magnitudes by an amount (AmF8i4w ~ 0.25) 
comparable to the shift between the observed LFs. The resulting internal distribution also 
has a larger width {a ~ 0.4) than the intrinsic one, due to the larger photometric errors. 

Our results show that the detected change in the C star brightness over NGC 4449 
field can be largely attributed to differences in completeness between the center and the 
most external regions. Accounting for this effect, the average magnitude of the C stars is 
constant within the errors (~ 0.05 mags). From the BD05 relation □, a change in magnitude 
of 0.05 corresponds to A[Fe/H] ^ 0.2. We interpret this as an upper limit to the metallicity 
variation over the field of view. This is consistent with st udies of metallicity gradients in 



other magellanic irregulars. For example, I Cole et al.l (120041 ) derive a metallicity gradient of 



0.05 dex X kpc ^ i n the LMC by c ompari ng the abundance s in th e inner disk and in 



the outer disk/spheroid (lOlszewskil 119931 ) . while iGrocholski et al.l (120061 ) find no metallicity 



gradients from spectroscopic studies of cluster stars in the LMC. 



4.2. RGB stars 

In the CMD of Fig. El the morphology of the RGB, at mF8i4w ^ 24 and mF555w ~ mF8i4W 
is connected to the properties of the stellar content older than ~ 1 Gyr. Despite the poorer 



^ BD05 use the Johnson-Cousins I-band. This is very similar to mF8i4W ( Sirianni et al.ll2005l ). and the 
small difference can be ignored for the purpose of the differential argument presented here. 
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time resolution with increasing look-back time, and the well known age-metallicity degener- 
acy, some constraints on the properties of the old stellar population can be inferred from an 
analysis of the RGB morphology. 

We derived the average RGB color as a function of mF8i4w by selecting stars with 
niF8i4W ^ 24, and then performing a Gaussian fit to the mpsssw — niF8i4w color distribution 
at different magnitude bins. The peak of the Gaussian fit is as red as mF555w ~ niF8i4W ~ 1-7 
at the RGB tip, and it is mF555w ~ mF8i4W ~ 1-45 at mF8i4w = 25, one magnitude below 
the tip. As expected, the RGB is significantly redder than in more metal-poor sta r-forming 



galaxies that we have previously studied with HST/ACS (lAloisi et al.ll2005l . 120071 ). 



In order to reveal the presence of age/metallicity gradients in NGC 4449, we performed 
a spatial analysis of the RGB morphology, following the same procedure as in Section 4.1. 
For each of the 28 rectangular regions identified in Fig. [71 we performed a Gaussian fit to the 
RGB mF555w — niF8i4w color distribution for bins of AmF8i4w = 0.25. Then we averaged for 
each region the colors derived in the four brightest bins, at 24 < mFsssw — niF8i4w < 25, i.e. 
one magnitude below the TRGB. The results of our analysis are presented in Fig. [T21 We 
plot along the ordinate the difference between the RGB color in each region and the average 
RGB color in the total field of view of NGC 4449; along the abscissa is the X coordinate in 
pixels. From top to bottom, the panels refer to regions of decreasing Y coordinate. Fig. [T2] 
shows that the RGB is bluer in the center than in the periphery of NGC 4449, with variations 
up to ~ 0.3 mag in the mF555w ~ mF8i4W color. 

This effect is also shown in the top panel of Fig. [131 where we plotted the observed 
color distributions of stars with 24 < mF8i4 < 25, for an external region ((1,1:4), (1:7,4), 
(7,1:4)), and for an internal region (3:5,2) of NGC 4449. The peaks of the Gaussian fits to 
the external and internal distributions differ by ~ 0.26 mag, and their errors are very small 
{a/y/N PS 0.001 mag) due to the large number of stars in each distribution 20,000). 
The Gaussian fit is broader for the internal region {a ~ 0.4) than for the external region 
(cr ^ 0.2). Also, the central region has a much broader tail of stars towards blue colors, due 
to the contamination from younger blue-loop and MS stars. The contribution from the MS 
+ blue-loop stars at the hot edge of the core He-burning phase is recognizable as a bump at 
V-I ^ 0. 

As was done for the C stars in Section 4.1, we performed Monte Carlo simulations to 
understand if the observed difference in the color distributions is intrinsic, or if it can be 
attributed to the larger crowding in the central regions of NGC 4449. We assumed an initial 
distribution, and drew Monte Carlo extractions from it with application of photometric errors 
and incompleteness. As a first guess for the initial distribution, we adopted a Gaussian with 
parameters mF555w ~ niF8i4w = 1.65 and a = 0.15. This has the same mean but somewhat 
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smaller a than the observed color distribution in the most external regions. The simulated 
color distributions for the external and internal regions were generated by applying the 
photometric errors and the completeness levels derived from artificial star experiments in the 
considered regions of NGC 4449. The results are presented in the central panel of Fig. [T31 
As observed, the width of the simulated distribution in the internal region is larger than in 
the external one because of the higher photometric errors. While the peak of the external 
simulated distribution is the same as that of the initial distribution, the peak of the internal 
distribution is blueshifted by an amount of ~ 0.06 mags. However, this is much less than the 
observed shift of ~ 0.3 mags. The simulations therefore show that the completeness and the 
photometric error variations over NGC 4449 field of view can account only in part 20 % 
) for the shift between the observed internal and external distributions. Thus the observed 
shift must be mainly due to an intrinsic variation of the stellar population properties. 

As a test, we performed a new simulation starting from a Gaussian with a peak as 
blue as mpsssw — niF8i4 = 1.45 and with a = 0.15. The result for the internal region is 
shown in the bottom panel of Fig. [131 Here we plot also the observed color distribution, 
for comparison. Both the peaks and the widths of the simulated and observed distributions 
are in good agreement. The only discrepancy is observed at the bluest color, where of 
course we do not reproduce the tail of the MS and post-MS stars. Our simulations suggest 
that the peak of the color distribution toward the center is intrinsically bluer than in the 
periphery. Assuming that the bluer peak in the center is due to a bluer RGB, the possible 
interpretations are 1) younger ages; 2) lower metallicities; or 3) lower reddening. If differential 
reddening is present within NGC 4449, we expect the central regions to be more affected 
by dust extinction than the periphery, and this would cause an even redder RGB in the 
center. Lower metallicities in NGC 4449 center are also very unlikely, since abundance 
determinations in galaxies show that metallicity tends to decrease from the center outwards 
or to remain flat (see Section 4.1). Thus, a bluer RGB would most likely indicate a younger 
population in the center of NGC 4449. Alternatively, the bluer peak observed in the center 
could be due to "contamination" by intermediate-age blue loop stars at the red edge of their 
evolutionary phase. But this too would imply the presence of a younger stellar population in 
the center than in the periphery of NGC 4449. More quantitative results on age /metallicity 
in NGC 4449 gradients will be derived through fitting of synthetic CMDs, and presented in 
a forthcoming paper (Annibali et al. 2008, in preparation). 
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A new distance determination 



The magnitude of the TRGB can be used to determine the distance of NGC 4449. The 
top panel of Fig. [H] shows the I-band LF of those stars in our final catalog that have V— I 
in the range 1.0-2.0. Here V and I are Johnson- Cousins magni tudes, obtained fro m our 



magnitudes in the ACS filter system using the transformations of lSirianni et al. (I20051) and 
applying a foreground extinction correction of E(B — V) = 0.019 ( Schlegel et al. 19981 ). The 
TRGB is visually identifiable as the steep increase towards fainter magnitudes at I ~ 24. At 
this magnitude, RGB stars start to contribute with a LF that increases roughly as a power 
law towards faint magnitudes. By contrast, the stars in the LF at brighter magnitudes are 
exclusively red supergiants and AGB stars. The drop at magnitudes fainter than I ^ 25.5 is 
due to incompleteness. 

To determine the TRGB magnitude we used the software and meth odology developed 
by one of us (R.P.v.d.M) and described in detail in ICioni et al.l (120001 ). A discontinuity 
produces a peak in all of the higher-order derivatives of the LF. We use a so-called Savitzky- 
Golay filter on the binned LF to obtain the first- and second-order derivatives. These are 
shown in the middle and bottom panel of Fig. [TH respectively. Peaks are indeed visible 
at the expected position of the TRGB. We fit these with Gaussians and find that the first 
derivative has a peak at Ii = 24.09, while the second derivative has a peak at I2 = 23.88. The 
reason for the difference between these magnitudes is the presence of photometric errors and 
bi nning in the analys is. This smooths out the underlying discontinuity. As shown in Fig. A.l 
of ICioni et al.l (I2OOOI ). this causes the first derivative to overestimate the magnitude of the 
TRGB, and the second derivative to underestimate th e magnitud e of th e TRGB. These 
biases can be explicitly corrected for as in Fig. A. 2 of ICioni et al.l (I2OOOI ). using a simple 
model for the true underlying LF and the measured width of the Gaussian peaks in the first- 
and second-order derivatives of the LF. After application of this correction we obtain the 
final estimates Itrgb,i = 23.99 and Itrgb,2 = 24.00 for the underlying TRGB magnitude, 
based on the first and second derivatives, respectively. The good agreement between these 
inde pendent estimates s hows that the systematic errors in the method are small, in agreement 
with ICioni et al.l (I2OOOI ) who adopted a systematic error AItrgb = ±0.02. The additional 
systematic error introduce d by the uncertaintie s in photometric zeropoints, transformations, 
and aperture corrections (jSirianni et al.l l2005l ) is AItrgb = ±0.03. The random error on 
AItrgb is very small due to the large number of stars detected in NGC 4449. It can be 
estimated using bootstrap techniques to be AItrgb = ±0.01. 

Our estimate of the TRGB magnitude, Itrgb = 24.00 ± 0.01 (random) ±0.04 (system- 
atic), can be compared to the abs olute magnitude of the TRGB, which was calibrated 
as a function of metallicity by, e.g., iBellazzini et al.l (120041 ). Adopting [M/H] = —0.52 for 
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NGC 4449 (based on the oxygen abundance given in Section 1) their cahbration (top panel 
of their Fig. 5) predicts Mjtrgb = —3.91. Comparison of different studies suggests that 
the systematic uncertainty in this prediction is ~ 0.15. This takes into account also the 
possibility that the RGB star metallicity is actually lower than that of the HIT regions 
(see discussion in Section 6). The implied distance modulus for NGC 4449 is therefore 
(m — M)o = 27.91 ±0.15, where we have added all sources of uncertainty in quadrature. 
This corresponds to D = 3.82 ± 0.27 Mpc. 

An alternative method for estimating the distance of NGC 4449 is through the average 
I-band magnitude of the carbon stars. Fig. [15] shows the I-band luminosity function of 
those stars in our final catalog that have V— I in the range 2.2-3.0. There is a well-defined 
peak, due to the horizontal "finger" of carbon stars seen in the CMD of Fig. O A Gaussian 
fit to the peak yields Icarbon = 23.59 ± 0.01, corrected for foreground extinction. Adopting 
[Fe/H] = —0.52 for NGC 4449 (based on the oxygen abundance given in Section 1, and 
assuming for simplicity that oxygen traces the iron content), the BD05 calibration in ([3]) 
predicts Mj^carbon = —4.48. The systematic uncertainty in this prediction is difficult to 
quantify. This is because carbon-star magnitudes are not well understood on the basis of 
stellar evolution theory (by contrast to the TRGB), because the dependence on stellar age 
or star formation history is poorly quantified, and because only a few empirical studies exist. 
We therefore adopt an uncertainty of ~ 0.2 mag, consistent with the discussion of Section 4.1 
The implied distance modulus for NGC 4449 is then (m — M)o = 28.07 ± 0.20, where we have 
added all sources of uncertainty in quadrature. This corresponds to D = 4.11 ± 0.38 Mpc. 
This is in agreement with the TRGB result, given the uncertainties. 

We can compare our result o f D = 3.82 ± 0.27 Mpc to the previous estimate of NGC 4449 



distance by iKarachentsev et al.l (120031 ). who inferred D = 4.2 ± 0.5 Mpc. This result was 
derived by applying the TRGB-luminosity method to HST/WFPC2 data. The extinction- 
corrected magnitude at which they detected the TRGB is Itrgb = 24.07 ± 0.26, which is 
consistent with our somewhat lower value of Itrgb = 24.00 ± 0.01 (random) ±0.04 (system- 
atic). Their somewhat larger distance is also due to the different value adopted for the 
absolute magnitude of the TRGB. They adopt Mj trgb = —4.05 (as appropriate for met al- 



poor systems), while we adopt Mi trgb = —3.91 from the more recent iBellazzini et al.l (120041 ) 



calibration (which takes into account the metallicity dependence of the TRGB luminosity). 



6. Comparison with models 



For a direct interpretation of the CMDs in terms of the stellar evolutionary phases, 
we have superimposed stellar evolutionary tracks for different metallicities on the mpsssw? 
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mF435W — niF555w, and mF8i4W, nipsssw — niF8i4W CMDs (Figs. [16] and [T7]). The tracks at 
Z=0.008, Z=0.004, and Z=0.0004 are the Padua stellar evolutionary tr acks (Fagotto et al. 



1994a , 1994b) transformed into the ACS Vegamag system by applying t helOriglia &: Leitherer 



((20001) code, and corrected for Galactic extinction (E(B - V) = 0.019. ISchlegel et al.l fll998f )) 



and distance modulus ((m — M)o = 27.91, see Section 5). T he Z=0.001 tracks wer e obtained 



from the Padua tracks through interpolation in metallicity (lAngeretti et al.l 120071 ). 



The metallicity range covered by the plotted tracks can account for the different pop- 
ulations potentially present in NGC 4449. The Z=0.008 and Z=0.004 tracks are suitable 
to account for the more metal rich population, given that the abundances derived in the 
nucleus and disk of NGC 4449 are almost one third of the solar value (see Section 1). The 
plotted tracks are for masses in the range 0.9-40 Mq. Lower mass stars would not have 
the time to reach visible phases within a Hubble time at the distance of NGC4449 (in these 
sets a 0.8 Mq star reaches the TRGB in 19 Gyr). The tracks are divided into three groups, 
namely low-mass stars (M < Mhgf), intermediate-mass stars (Mhcf < M < Mup) and high- 
mass stars (M > Mup). The subdivision is made according to the critical mass at which the 
ignition of the central fuel (either helium or carbon) starts quietly depending on the level 
of core electron-degeneracy. In the adopted tracks, the value of Mhgf depends shghtly on 
metalhcity, being equal to 1.7, 1.8 and 1.9 for Z=0.0004, Z=0.004 and Z=0.008, respec- 
tively. The value of Mup is between 5 and 6 Mq. For low- mass stars, we have displayed in 
Figs. [16] and [17] only the phases up to the TRGB in order to avoid excessive confusion. 

At the distance of NGC 4449, the MS is sampled only for stars with M > 3Mq, and 
corresponds to the vertical lines whose mF435w — niFsssw, niFsssw — niF8i4w colors go from 
~ —0.3 to ~ from the 40 Mq to the 3 Mq track. The turnoff is recognizable as a small blue 
hook on each evolutionary track of Fig. [161 The almost horizontal blue loops correspond 
to the later core helium-burning phase for intermediate- and high- mass stars, while the 
bright vertical lines in the red portion of the CMD describe the AGB phase. For low- 
mass stars, our data sample only the brighter red sequences of the RGB, all terminating 
at the TRGB with approximately the same luminosity. The models considered here do not 
include the horizontal feature of the carbon stars in the TP-AGB phase, recognizable in the 



H1F814W5 mpsssw — niF8i4W CMD at mF555w — niF8i4W ^ 2, but see iMarigo fc Girardi] (120071 ) 
for new appropriate models. 

The comparison of the CMDs with the tracks in Fig. [16] shows that the blue plume is 
populated by high- and intermediate-mass stars on the MS, and high-mass stars at the hot 
edge of the core helium-burning phase. The red plume samples bright red SGs and AGB 
stars. The faint red sequence at mF8i4w ^ 24, featured in Fig. [TT] is due to low-mass stars 
in the RGB phase. 
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The presence of high-, intermediate- and low-mass stars in the CMD testifies that young, 
intermediate-age and old stars (several Gyrs) are present at the same time in NGC 4449. 
In particular, the fact that we sample stars as massive as 40 Mq implies that the SF was 
active 5 Myr ago. Furthermore, the absence of significant gaps in the CMD suggests that the 
SF has been mostly continuous over the last 1 Gyr. The time resolution gets significantly 
poorer for higher look-back times because of both the intrinsic degeneracy of the tracks and 
the large photometric errors at faint magnitudes; thus small interruptions in the SF can be 
easily hidden in the CMD for ages older than > 1 Gyr. 

The metallicity [M/H] = —0.52 derived i n Section 1 for NGC 44 49 corresponds to a 



metal fraction Z=0.005, adopting Z© = 0.017 (iGrevesse fc Sauvallll998l ). However this value 



refers to abundance determinations in HII regions, thus it is likely to refiect the metallicity 
of the youngest stars. Fig. [12] and [T7I show that the Z=0.004 tracks are in good agreement 
with all the phases of the empirical CMDs. The Z=0.0004 and Z=0.001 tracks are too blue 
to account for the observed RGB feature. So, if there is no significant extinction intrinsic to 
NGC 4449, then stars with such low metallicities do not account for a significant fraction of 
the stellar population in NGC 4449. Allowing for an age spread from 1 Gyr to a Hubble time, 
the Z =0.004 tracks are in good agreement with both the color and the width of the observed 



RGB. iHill et all fll998h derived an internal reddening E(B - V) ^ 0.18 from the Ha/R/S 
ratio measured in NGC 4449 HII regions. Adopting this value, the RGB is consistent with 
metallicities as low as Z=0.001. However, this extinction value is appropriate for the young 
star forming regions, and we expec t a lower e xtinct ion in the regions where older stars reside 
(as demonstrated for the LMC by IZaritsky Jl999h l Either way, the Z=0.008 models seem 



slightly too red, especially in the AGB phase. Although there are some uncertainties in the 
models of this phase, we do believe that this discrepancy indicates that metallicities higher 
than Z=0.008 are ruled out for this galaxy. 

By superimposing the stellar tracks on the empirical CMDs, we can attempt a selec- 
tion of the photometrized stars according to their age. In Fig. [18] we show the four regions 
selected on the mF8i4W5 mF555W — niF8i4w CMD that correspond to different stellar masses, 
namely M> 20 M©, SM© < M < 20Mq, I.8M0 < M < SMq, and M < I.8M0. These re- 
gions roughly define the loci of very young stars, with ages < 10 Myr, young stars, with 
10 Myr < age < 100 Myr, intermediate- age stars, with ages between ~ 100 Myr and 1 Gyr, 
and old stars, with ages > 1 Gyr. The spatial distribution of the four groups of stars is shown 
in Fig. [191 We notice that old stars are homogeneously distributed over the galaxy, except 
for the central regions, where the high crowding level makes their detection more difficult. 
As we approach younger ages, the distribution becomes more and more concentrated. Stars 
with ages between 10 and 100 Myr are highly clustered in an S shaped structure centered 
on the galaxy nucleus. This could be a bar, which is a common feature among Magellanic 
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irregular galaxies. Stars with ages between 10 and 100 Myr also clearly outline the symmet- 
ric structure around the resolved cluster-like object (see Fig. [7]) that was already discussed 
in Section 4. The same structure is visible to a lesser extent in the spatial distribution of 
stars with ages 100 Myr - 1 Gyr. The resolved cluster-like object itself has a clear RGB (not 
shown here), and it is therefore not a young super star cluster. Very young stars, with ages 
younger than 10 Myr, are strongly clustered and detected only in the very central regions of 
NGC 4449. Such young stars should be able to ionize the surrounding interstellar medium 
and produce HII regions. Indeed, Fig. [201 where we plotted the F658N (Ha) image together 
with the positions of the stars younger than 10 Myr, shows that this is the case, with a strong 
correlation between the position of the stars younger than 10 Myr and the HII regions. This 
indicates that most of the emission is due to photoionization rather than to shocks, due to, 
e.g., supernovae explosions. 



7. Conclusions 

We have acquired HST/ACS imaging in the F435W (B), F555W (V), F814W (I) and 
F658N (Ha) filters of the Magellanic starburst galaxy NGC 4449 in order to infer its star 
formation history and understand the properties of the observed global starburst. In this 
paper we present the B, V and I photometry of the resolved stars. We detect 299,014 objects 
in the (B,V) CMD, 402,045 objects in the (V,I) CMD, and 213,099 objects with a measured 
magnitude in all the three bands. The derived CMDs span a magnitude range of ~ 10 
mag, and sample both the young and the old resolved stellar population in NGC 4449. We 
also detected several candidate clusters (at least ~ 40, some of which look like very massive 
globular clusters) and background galaxies in our images. 

We derived a new distance from the TRGB method. The TRGB is detected at a 
Johnson-Cousins I magnitude of Itrgb = 24.00 ± 0.04. At the metallicity of NGC 4449, the 
TRGB is expected at an absolute magnitude of Mj trgb = —3.91, with a systematic error of 
~ 0.15 mag. This provides a distance modulus of (m — M)o = 27.91 ± 0.15, i.e. a distance 
of 3.82 ± 0.27 Mpc. We used also the alternative method of the carbon-star luminosity, and 
found a distance of 4.11 ± 0.38 Mpc, which is consistent with the result from the TRGB 
method. Our distance determinations are consistent within the errors with the value of 
D= 4.2 ± 0.5 Mpc previously provided by Karachentsev et al. (2003). 

In the CMDs of NGC 4449 we observe a well defined blue plume (MS and post-MS 
stars) and red plume (red SGs and AGB stars), the horizontal tail of the carbon stars in 
the TP- AGB phase (in the I, V— I CMD), and a prominent RGB. The presence of all these 
evolutionary features implies ages up to at least 1 Gyr and possibly as old as a Hubble time. 
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The comparison of the observed CMDs with the Padua stellar evolutionary tracks, corrected 
for the derived distance modulus and foreground extinction, shows that stars as massive 
as 40 Mq are present in NGC 4449. Such high masses imply that the star formation was 
active 5 Myr ago, and possibly it is still ongoing. The absence of significant gaps in the 
CMDs suggests also that the star formation has been mostly continuous over the last 1 Gyr. 
However, interruptions in the star formation can be easily hidden in the CMD for ages older 
than > 1 Gyr, because of the intrinsic degeneracy of the tracks and the large photometric 
errors at faint magnitudes. The presence of a prominent RGB testifies that NGC 4449 hosts a 
population possibly as old as several Gyrs or more. However, the color-age degeneracy of the 
tracks with increasing look-back time (at a given metallicity, large age differences correspond 
to small color variation in the RGB), and the well known age-metallicity degeneracy, prevent 
us from establishing an exact age for the galaxy. We will derive a better age estimate with 
the synthetic CMD method, when we'll study the detailed SFH (Annibali et al. 2008, in 
preparation). 



Abundance es timates in NGC 4449 HII regions provide 12 + log(0/H)f« 8.31 (ITalent 



1980l : lMartinlll997l ). whi ch corresponds to [O/H] = —0.52 if we assume a solar abundance of 



12 -|- log(O/H)0 = 8.83 (iGrevesse &: Sauvallll998l ). These measures are biased toward regions 



where the interstellar medium has been significantly reprocessed, and thus are likely to reflect 
the metallicity of the youngest generation of stars. We expect lower metallicity for the oldest 
stars in NGC 4449. Interestingly, though, the Z=0.004 stellar evolutionary tracks (the closest 
in the Padua set t o the metallicity of NGC 4449), corresponding to log(Z/Z0) = —0.63 if we 



adopt Zq = 0.017 (IGrevesse &: Sauvallll998l ). seem to be in very good agreement with all the 
features observed in the empirical CMDs, if we assume that there is not significant extinction 
intrinsic to NGC 4449. In particular, the next lower metallicity tracks (at Z = 0.001) are 
definitively too blue to account for the observed RGB colors, implying that the bulk of the 
stellar population older than 1 Gyr was already enriched in metals. 

We investigated the presence of age and metallicity gradients in NGC 4449. To this 
purpose, we divided the total galaxy's field of view into 28 rectangular regions of ~ 1 kpc^ 
area, and derived the CMDs for the different regions. The CMD morphology presents a 
significant spatial dependence: while the RGB is detected over the whole field of view of the 
galaxy, the blue plume, red plume and blue-loop stars are present only in the more central 
regions, indicating that the stellar population is younger in the center than in the periphery 
of NGC 4449. 

We also studied the spatial behavior of the carbon-star luminosity and of the RGB 
color. Once the effect of incompleteness and photometric errors is taken into account, the 
average magnitude of the C stars turns out to be constant within the errors (~ 0.05 mags). 
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This gives an upper limit of ~ 0.2 dex on the metalhcity variation over the field of view of 
NGC 4449. On the other hand, we find that the RGB is intrinsically bluer in the center 
than in the periphery of the galaxy. Bluer RGB colors can be due to younger and/or more 
metal poor stellar populations. However, as spectroscopic-based abundance determinations 
in galaxies show that metallicity tends to decrease from the center outwards, or to remain 
constant, we interpret this as the result of a younger, and not more metal poor, stellar 
population in the center of NGC 4449. 

With the help of the Padua tracks, we identified in the observed mF8i4w, nipsssw — niF8i4w 
CMD four different zones corresponding to different stellar masses, namely M> 2OM0, 
5M0 < M < 2OM0, I.8M0 < M < 5M0, and M < I.8M0. These regions roughly define the 
loci of stars with age < 10 Myr, 10 Myr < age < 100 Myr, 100 Myr < age < 1 Gyr, and 
age > 1 Gyr. Low-mass old stars are homogeneously distributed over the galaxy's field of 
view, with the exception of the central regions, where the high crowding level makes their 
detection more difficult. As we approach younger ages, the spatial distribution of the stars 
becomes more and more clustered. Intermediate- age stars (100 Myr < age < 1 Gyr) are 
mostly found within ^ 1 kpc from the center. Stars with ages between 100 Myr and 10 Myr 
are found in an S-shaped structure centered on the galaxy nucleus, and extending in the 
North-South direction up to 1 kpc away from the center. This could be a bar, which is a 
common feature among Magellanic irregular galaxies. Stars younger than 10 Myr are very 
rare, and found only in the galaxy nucleus and in the North arm of the S shape structure. 
The comparison with the Ha image shows a tight correlation between the position of the 
stars younger than 10 Myr and the HH regions, indicating that we have identified the very 
massive and luminous stars that ionize the surrounding interstellar medium. 

One of the many star clusters visible in our image is of particular interest. This cluster 
on the West side of the galaxy is surrounded by a symmetric structure that is particularly 
well outlined by stars with ages in the range of 10-100 Myr (see Fig. fT^ . This structure could 
be due to tidal tails or spiral-like feature associated with a dwarf galaxy that is currently 
being disrupted by NGC 4449. The cluster could be the remnant nucleus of this galaxy. It 
is resolved into red stars, and it has a significant ellipticity (see Fig. [7]). This is reminiscent 
of the star cluster u Cen in our own Milky Way, which has also beeri suggested to be the 
remnant nucleus of a disrupted dwarf galaxy (iFreeman Sz Rodgers Il975h . More details about 



this object will be presented in a forthcoming paper (Aloisi et al. 2008, in preparation). 

Quantitative information on the star formation history of NGC 4449 is fundamental 
in order to understand the connection between the global starburst observed and processes 
such as merging, accretion and interaction. The star formation history of NGC 4449 will be 
derived in a forthcoming paper through the synthetic CMD method, which is based on stellar 
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evolutionary tracks, and is able to fully account for the effect of observational uncertainties, 
such as photometric errors, blending and incompleteness of the observations. 

Support for proposal #10585 was provided by NASA through a grant from STScI, which 
is operated by AURA, Inc., under NASA contract NAS 5-26555. We thank Livia Origlia for 
providing the photometric conversion tables to the ACS Vegamag system. 
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Fig. 1.— Four-color (F435W (B), F555W (V), F814W (I) and F658N (Ha)) composite image 
of NGC 4449 showing a field of view of ~ 380 x 200 arcsec^, obtained from the mosaic of 
the two different ACS pointings, the B, V, I and Ha data are shown in blue, green, red and 
pink, respectively. The image orientation is also indicated. 
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Fig. 2.— Luminosity functions (LFs) in the F435W, F555W, and F814W filters. Tlie solid 
line is the LF obtained if the star detection is performed independently on the three different 
images. The dot-dashed line refers to a deeper photometry, obtained by detecting the stars 
in the sum of the F435W, F555W, and F814W images. 
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Fig. 3. — Distribution of the DAOPHOT parameters a, sharpness and ^is a function of 
the magnitude. The a distributions for the F435W (B), F555W (V) and F814W (I) fihers 
are shown in the three upper panels. The sharpness and distributions are shown only for 
the F814W filter, but look very similar in F435W and F555W. 
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Fig. 4. — Input minus output magnitude versus input magnitude of the artificial stars in the 
F435W (top), F555W (middle) and F814W (bottom) filters. For each filter, we simulated 
half of a million stars. The lines superimposed on the diagrams represent the mean Am and 
the ±1 standard deviations. 
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Fig. 5. — niF555w, inF435w — inF555w CMD of the 299,014 stars matched between the B and 
V catalogs. The 90 % (sohd hne) and 50 % (dashed hne) completeness levels are indicated 
on the CMD. The average size of the photometric errors, as derived from the artificial star 
experiments, is indicated at different mpsssw magnitudes. The color error is evaluated at 

I11F435W — niF555W = 1- 
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Fig. 6. — niF8i4w, mF555w — niF8i4w CMD of the 402,045 stars matched between the V and 
I catalogs. The 90 % (sohd hne) and 50 % (dashed hne) completeness levels are indicated 
on the CMD. The average size of the photometric errors, as derived from the artificial star 
experiments, is indicated at different mF8i4w magnitudes. The color error is evaluated at 

I11F555W — niF814W = 1- 
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Fig. 7. — Mosaicked image of NGC 4449 where the dotted hnes indicate the adopted subdi- 
vision into (4 X 7) regions. The insets show blowups of ~ 5 x 5 arcsec^ regions to provide 
an idea of both the amount of crowding and the resolution of our data. At the distance of 
NGC 4449, each region samples an area of ~ 1 kpc^. 
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Fig. 8. — niF555w versus mF435w — nipsssw CMDs for the 28 regions selected in the field of 
NGC 4449. The CMDs are represented through Hess diagrams. The number density of stars 
in the mpsssw? mF435W ~ mpsssw plane increases from blue to red. We indicate the average 
photometric errors as derived from the artificial star experiments for each region. The color 
error is evaluated at mF435w — nipsssw = 1- The 90 % (solid line) and 50 % (dashed line) 
completeness levels are indicated on the CMDs of the central column. 
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Fig. 9. — niF8i4w versus mpsssw — niF8i4w CMDs for the 28 regions selected in the field of 
NGC4449. The CMDs are represented through Hess diagrams. The number density of stars 
in the mF8i4W) mpsssw ~ mF8i4W plane increases from blue to red. We indicate the average 
photometric errors as derived from the artificial star experiments for each region. The color 
error is evaluated at mpsssw — niF8i4w = 1- The 90 % (solid line) and 50 % (dashed line) 
completeness levels are indicated on the CMDs of the central column. 
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Fig. 10. — Variation of the carbon star luminosity as a function of position in the field 
of NGC 4449. The quantity Am/^car&on along the ordinate is the difference between the 
carbon star F814W magnitude measured in a specific region, and the carbon star magnitude 
averaged over the whole field of view of NGC 4449. Along the abscissa is the X coordinate 
in pixels. From top to bottom, the panels refer to regions of decreasing Y coordinate. 
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Fig. 11. — Top panel: observed C star luminosity functions (LFs) for an external region 
((1,1:4), (1:7,4), (7,1:4)) of NGC 4449 (shaded histogram), and for the most internal one 
(4,3). The peaks of the external and internal distributions, as derived from Gaussian fits, 
are mF8i4 = 23.64 and mF8i4 = 23.37, respectively. The external distribution is narrower 
(cf814W ~ 0.3) than the internal one (o"f814W ~ 0.5). Bottom panel: simulated distribution 
for the C stars. The Gaussian curve is the assumed initial mF8i4 distribution for the C stars 
(mF8i4,o = 23.64, (Jfsmw = 0.3). The shaded and thick line histograms represent the resulting 
distributions in the external and internal regions, after errors and incompleteness have been 
applied. The peak of the simulated internal distribution has a shift of AmF8i4w ~ —0.25 
with respect to the peak of the Gaussian curve. Its width (a ^ 0.4) is larger than the 
Gaussian, due to the larger photometric errors in the center of NGC 4449. 
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Fig. 12. — mF555w ~ niF8i4w color variation of the RGB as a function of position in NGC 4449 
field. The RGB colors were computed performing an average between 24 < mF8i4w ^ 25, i.e. 
within one magnitude below the TRGB. The variation AmF555w ~ niF8i4w was computed 
with respect to the RGB color in the total field of view of NGC 4449. Along the abscissa is 
the X coordinate in pixels. Prom top to bottom, the panels refer to regions of decreasing Y 
coordinate. 
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Fig. 13. — Top panel: Observed color distributions of stars with 24 < mF8i4 < 25 
mag, for an external region ((1,1:4), (1:7,4), (7,1:4), shaded histogram), and for an in- 
ternal region ((3:5,2) thick line) of NGC 4449. The distributions, which peak at red 
(mF555w — mF8i4w > 1) colors, denote the presence of old low-mass stars in the RGB phase. 
The dotted lines indicate the position of the peaks (as derived from a Gaussian fit) of the 
external and internal distributions, at mF555w ~ iiif814 = 1-65 and mF555w ~ niF8i4 = 1.39, 
respectively . Central panel: simulated color distributions. The Gaussian curve is an 
assumed intrinsic distribution, with peak and standard deviation as indicated in the panel. 
The shaded and thick line histograms are the simulated distributions in the external and in- 
ternal regions, respectively. They were generated through Monte Carlo extractions from the 
assumed Gaussian, and applying photometric errors and completeness levels as derived from 
artificial star experiments in the considered regions. The dotted lines indicate the peaks of 
the simulated distributions. While the peak of the external distribution is the same as that 
of the initial Gaussian, the peak of the internal distribution is blueshifted by ^ 0.06 mags. 
Bottom panel: The thick line is the observed color distribution for the internal region 
(3:5,2). The shaded histogram is a new simulated distribution for the central region drawn 
from a Gaussian with a bluer peak. The peaks of the observed and simulated distributions, 
as indicated by the dotted line, are now the same. 
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Fig. 14. — Top panel: extinction corrected /-band luminosity function (LF) of stars with V— I 
in the range 1.0-2.0 (sohd). The solid line is our final estimate of the TRGB magnitude, 
obtained as described in the text. Middle panel: first-order derivative of the LF (solid), 
with the best-fitting Gaussian overplotted (dashed). The vertical line indicates the center 
of the Gaussian, which provides a slightly biased overestimate of the TRGB magnitude. 
Bottom panel: second-order derivative of the LF (solid), with the best-fitting Gaussian 
overplotted (dashed). The vertical line indicates the center of the Gaussian, which provides 
a slightly biased underestimate of the TRGB magnitude. The vertical scales in all panels 
are normalized arbitrarily. 
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Fig. 15. — Extinction corrected J-band luminosity function (LF) of stars with V— I in the 
range 2.2-3.0 (sohd), with the best-fitting Gaussian overplotted (dashed). The vertical line 
indicates the center of the Gaussian, which is an estimate of the average /-band magnitude 
of the carbon stars in NGC 4449. The vertical scale is normalized arbitrarily. 
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Fig. 16.— mpsssw, mF435W - nipsssw CMDs (left panels) and nipsuw, mpsssw - hifsmw 
CMDs (right panels) of NGC 4449, with superimposed stellar evolutionary tracks for different 
metallicities, as indicated in the panels. The CMDs are represented through Hess diagrams, 
with the number density of stars increasing from blue to red. Points have been used to 
plot the brightest stars. The Z=0.008, Z=0.004 and Z=0.0004 models are the Padua stellar 
evolutionary tracks (Fagotto et al. 1994a, 1994b) converted into the ACS filters by applying 
bolometric corrections, and adopting a foreground E(B — V) = 0.019 and distance modulus 
(m — M)o = 27.91. The Z= 0.001 tracks were obt ained from the Padua tracks through linear 
interpolation in metallicity (lAngeretti et al.ll2007l ). We displayed the following stellar masses: 
(left to right): 40, 30, 20, 12, 9, 7, 5, 4, 3, 2 and 1 Mq. The corresponding lifetimes (which 
slightly depend on metallicity) are: 5, 7, 10, 20, 35, 56, 110, 180, 370, 1120 and 8410 Myr, 
respectively. No extinction was assumed intrinsic to NGC 4449. 
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Fig. 17. — mF8i4w, niF555w — niF8i4w CMD of NGC 4449 showing the RGB feature, with 
superimposed the Padua low-mass stellar evolutionary tracks for different metallicities, as 
indicated in the panels. The CMDs are represented through Hess diagrams, with the number 
density of stars increasing from blue to red. The foreground extinction and the distance 
modulus adopted to display the tracks are E(B — V) = 0.019 and (m — M)o = 27.91. We 
plot the following stellar masses: (left to right): 1.8, 1.6, 1.4, 1.2, 1 and 0.9 Mq. The 
corresponding lifetimes (which slightly depend on metallicity) are: 1300, 1630, 2520, 4260, 
8410, 12500 Myr, respectively. No extinction was assumed intrinsic to NGC 4449. 
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Fig. 18. — mF8i4w, niF555w — niF8i4w CMD of NGC 4449 (402,045 stars), where we have 
identified the loci approximately corresponding to very young stars (age < 10 Myr), young 
stars (10 Myr < age < 100 Myr), intermediate-age stars (100 Myr < age < 1 Gyr), and old 
stars (age > 1 Gyr). On the same CMD, we plotted the stellar evolutionary tracks for the 
masses that indicatively define the age boundaries. 
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Fig. 19. — Spatial distribution of very young (age < 10 Myr), young (10 Myr < age < 100 
Myr), intermediate-age (100 Myr < age < 1 Gyr), and old stars (age > 1 Gyr) in NGC 4449. 
Stars were selected according to their position in the mF8i4w, niF555w ~ mF8i4W CMD, as 
shown in Figure [18l Notice that crowding prevents us from detecting old stars in the central 
regions. The red X marks the position of the resolved cluster-like object for which a blow-up 
is shown in Figure [71 This object is surrounded by a bi-symmetric structure that is clearly 
outlined by stars with ages from 10-100 Myr, and to a lesser extent by stars of ages 100 Myr 
- 1 Gyr. 




Fig. 20. — ACS-F658N (Ha) image of NGC 4449 plotted as a greyscale map. The positions 
of the stars younger than f» 10 Myr are superimposed using red dots. 



